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Abstract. In the periodic orbit quantization of physical systems, usually only the leading-order h contri- 
bution to the density of states is considered. Therefore, by construction, the eigenvalues following from 
semiclassical trace formulae generally agree with the exact quantum ones only to lowest order of h. In 
different theoretical work the trace formulae have been extended to higher orders of h. The problem re- 
mains, however, how to actually calculate eigenvalues from the extended trace formulae since, even with h 
corrections included, the periodic orbit sums still do not converge in the physical domain. For lowest-order 
semiclassical trace formulae the convergence problem can be elegantly, and universally, circumvented by 
application of the technique of harmonic inversion. In this paper we show how, for general scaling chaotic 
systems, also higher-order h corrections to the Gutzwiller formula can be included in the harmonic in- 
version scheme, and demonstrate that corrected semiclassical eigenvalues can be calculated despite the 
convergence problem. The method is applied to the open three-disk scattering system, as a prototype of a 
' chaotic system. 

O ' 

PACS. 03.65.Sq Semiclassical theories and applications 

m ' 
O 

CN 1 Introduction 

O . 

The relation between the eigenvalue spectrum of a quantum system and the periodic orbits of the corresponding 
classical system is a question of fundamental importance for both integrable and chaotic dynamical systems. The well- 
C ■ established Gutzwiller trace formula for classically chaotic systems and its analogue for integrable systems, 

J> \ the Berry- Tabor formula , provide the semiclassical density of states in terms of a sum over all periodic orbits of 
the system. However, each trace formula is only the leading-order term of an expansion of the exact density of states 
in powers of H, and therefore in general the resulting semiclassical eigenvalues are only approximations to the exact 
quantum ones. In recent years, two basic methods have been developed for determining higher-order h corrections 
to the Gutzwiller trace formula in terms of periodic orbit contributions [||[^,||||,|H],|ll]]. Unfortunately, even with h 
corrections included, the trace formulae usually suffer from being divergent in the region where the physical eigenvalues 
or resonances are located. For specific systems, higher-order H corrections to the semiclassical eigenvalues have explicitly 
been calculated by cycle expansion techniques [p|1s|,|lo||. However, this method is applicable only to systems with special 
features, namely to hyperbolic systems with a known complete symbolic dynamics. 

Recently, it has been demonstrated how the convergence problems of the semiclassical trace formulae can be 
circumvented by the application of harmonic inversion techniques [ |T2[|T3| ,|l4|]. The harmonic inversion method is capable 
of extracting semiclassical eigenvalues from a finite set of periodic orbits with very high precision and resolution. In 
contrast to other semiclassical methods, harmonic inversion does not require any special properties of the system, 
and can therefore be applied to a wide range of physical systems. In Refs. a general procedure has been 

developed for including higher-order h corrections to the trace formulae in the harmonic inversion scheme. So far, this 
method has only been tested for an integrable system, viz. the circle billiard. The general procedure, however, does 
not depend on the type of the underlying classical dynamics, and is applicable also to chaotic systems. In this paper, 
we demonstrate how the method works for chaotic systems, and apply it to the open three-disk scatterer, which has 
become a standard example for the semiclassical quantization of chaotic systems [ jj^[lj| ,[l7[[[8> 19, 2(i|,pl|. 



The harmonic inversion method is used in two directions: First, we carry out a harmonic analysis of the spectrum 
of the differences between the exact (complex) quantum eigenvalues and the semiclassical resonances of the three-disk 
system. We show how this enables one to determine, for each orbit, the first-order h correction term (and, in principle, 
all higher-order correction terms) to the Gutzwiller formula. We confirm our results by comparing with the values 
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calculated by a specialization of an analytical approach developed by Vattay and Rosenqvist |8|,|| for two-dimensional 
billiards [p_0| . Second, we take the analytical correction terms to the Gutzwiller formula and compute, using the 
classical periodic orbit data and harmonic inversion, the first-order h corrections to the semiclassical resonances of the 
three-disk system. Thus we illustrate that corrected semiclassical eigenvalues can be calculated by harmonic inversion 
despite the convergence problems of the trace formulae. We compare the zeroth and first-order approximations to the 
resonances with the exact quantum eigenvalues, and can quantitatively assess the increase in accuracy produced by 
including the next-order corrections. 



2 Higher-order H corrections to Gutzwiller's trace formula 

Gutzwiller's trace formula for chaotic systems gives a semiclassical approximation to the response function (i.e., the 
trace of the Green's function) of a quantum system in terms of the periodic orbits of the corresponding classical system. 
The semiclassical response function consists of a smooth background, and an oscillating part g(E) = g(E) + g osc (E), 
where the oscillating part is given by P,|22]] 



Sp 7T 



(1) 



The sum in ([!]) runs over all periodic orbits (po) of the system, including multiple traversals. Here, T po and S po 
are the period and the action of the orbit, M po and /i po denote the monodromy matrix and the Maslov index, and 
the repetition number r counts the traversals of the underlying primitive orbit ("primitive" means that there is no 
sub-period). The semiclassical density of states p(E) is related to the response function via 

p(E) = --lmg(E) . (2) 

In general, the semiclassical eigenvalues or resonances obtained from the Gutzwiller formula agree with the exact 
quantum ones only in leading order of h. To improve the accuracy of the semiclassical eigenvalues, higher-order h 
correction terms to the Gutzwiller formula have to be included. Two different methods for the calculation of such 
higher-order h terms have been derived for chaotic systems, one by Gaspard and Alonso and the other by Vattay 
and Rosenqvist [Bp, 10 . The latter method has been specialized to two-dimensional chaotic billiards in Ref. ]ic|] . An 



extension of the method of Gaspard and Alonso has recently been published in Ref. . We will adopt the method of 
Vattay and Rosenqvist to compute the first-order h corrections to the semiclassical resonances of the open three-disk 
scatterer. 

Vattay and Rosenqvist give a quantum generalization of the Gutzwiller formula, which is of the form 

g(E) = g(E) + I ]T]T (t p (E) i /^^ j f] (R l p (E)) r exp (^S P (E)\ . (3) 

p l \ / r=l ^ ' 

The first sum runs over all primitive periodic orbits; T p and S p are the traversal time and the action of the periodic 
orbit, respectively. The sum over r corresponds to multiple traversals of the primitive orbit. The quantities R l p are 
associated with the local eigenspectra determined by the local Schrodinger equation in the neighbourhood of the 
periodic orbits. An expansion of the quantities R l p in powers of h, 

R\E)^e X J±(^Yc^ 



v m=0 



m 



(O (l + ^Cf > + ...) , (4) 



« exp | 

yields the H expansion of the generalized trace formula (^). For two-dimensional hyperbolic systems, the zeroth-order 
terms are given by 



exp 



where p p and A p are the Maslov index and the expanding stability eigenvalue (i.e., the stability eigenvalue with an 
absolute value larger than one) of the orbit, respectively. By summation over /, the Gutzwiller trace formula is regained 



Kirsten Weibert et al.: Higher-order h corrections in the semiclassical quantization of chaotic billiards 



3 



as zeroth-order approximation to Eq. (||). If the zeroth-order terms do not depend on the energy, as is the case for 
billiard systems, the first-order correction to the Gutzwiller formula can be written as 

^)4^^^f eXp M ^ (1) expfe po (£)) , (6) 

po / ^ ' 

where the first sum in Eq. (Q) now runs over all periodic orbits, including multiple traversals, and r is the repetition 
number with respect to the underlying primitive orbit. 

An explicit recipe for the calculation of the correction terms Cj for two-dimensional chaotic billiards was given 
in Ref . |10| . The correction terms must in general be calculated numerically from the periodic orbit data. A numerical 
code which determines the first-order corrections for two-dimensional chaotic billiards can also be found in Ref. [fiof . 
We have used that code to compute the correction terms Cj for the three-disk system for a comparison with the 
correction terms calculated by harmonic inversion. 



3 The open three-disk scatterer 

As a model system for the calculation of higher-order h corrections to the Gutzwiller formula by harmonic inversion, 
we consider the open three-disk system, which consists of three equally spaced hard disks of unit radius. This system, 
in particular the case of the relatively large disk separation d — 6, has served as an archetype for the application of 
semiclassical quantization techniques in many investigations in r ecent years @grail9[|g,|l|. We will consider the 
case d = 6, as well as the small separation d = 2.5. In our calculations, we make use of the symmetry reduction of the 
three-disk system introduced in Refs. |18,2^] and concentrate on states of the A\ subspace. 

As for all billiar d systems, the shape of the periodic orbits in the three-disk system is independent of the wave 
number k = v2mE/H, and the action scales as 

S/h = ks, (7) 

where the scaled action s is equal to the physical length of the orbit. We consider the density of states as a function 
of the wave number 

p(k) = --Im g(k) , (8) 

with a scaled response function g(k). Since the wave number k is proportional to ftr 1 , it can be considered as an 
effective Planck constant, 

k = K-+ . (9) 
The h expansion of the exact quantum response function can therefore be written as a power series in k : 

g(k) = g(k)+g osc (k) (10) 

with 

oo oo 1 

g osc (k) = = j2 F E4o eiSpofc - (n) 

n— n— po 

The second sum runs over all periodic orbits including multiple traversals. The zeroth-order amplitudes Apo correspond 

to the Gutzwiller formula, whereas for n > 0, the amplitudes AplJ give the n th -order corrections g n (k) to the response 
function. 

Applying the Gutzwiller trace formula to the (symmetry reduced) three-disk system yields for the zeroth-order 
amplitudes in Eq. (0) (A\ subspace) 



42 = -iE 



-if Mpo 



po 



r|det(M po -l)| 1 /2 



where M po is the monodromy matrix of the orbit, l s is the symbol length, s po the scaled action, and A po the expanding 
stability eigenvalue of the orbit. The Maslov index /i po for this system is given by 21 s . The quantity r designates the 
repetition number with respect to the corresponding primitive orbit. The first-order amplitudes of the h expansion 
( p"l| ) following from Eq. (^|) read 

P ° " r ^ |A po |V2 A * o 2hk ■ [U) 



4 



Kirsten Weibert et al.: Higher-order h corrections in the semiclassical quantization of chaotic billiards 



Since the terms C\ are proportional to the momentum hk, as was shown in Ref. |10[], the amplitudes are independent 

of the scaling parameter k. The correction terms have to be determined numerically. We use the code developed 
by Rosenqvist and Vattay |T^,Q. The code requires the flight times between the bounces and the reflection angles 
as an input. These parameters have to be calculated numerically for each periodic orbit. As the contributions to the 
amplitude (|l3| ) for different I are proportional to |A po | _ ' - ^, the sum over / converges fast if the absolute value of the 
stability eigenvalue A po is large. For most orbits, the leading term I = turns out to be already sufficient. It is only 
for the very shortest orbits that terms of higher order in I have to be included to ensure convergence of the sum to 
within, say, 3 significant digits. 



4 Harmonic analysis of the quantum spectrum 

4.1 Theory 

In Refs. it was demonstrated that the amplitudes A p l) of the h expansion (|ll) can be obtained by a harmonic 

inversion analysis of the exact quantum spectrum. The general procedures do not depend on any special properties 
of the system and can be applied to both integrable and chaotic systems. In Ref. jl6|, they were tested for the circle 
billiard, as an example of an integrable system. We will now use the same procedures for the open three-disk system, 
as a representative of a chaotic system. 

We start by briefly recapitulating the main ideas of the procedures developed in Refs. Jl4] , [l6| . The exact quantum 
mechanical response function, in terms of the wave number k, can be written as 

*-<*>- ErngW (14) 

3 J 

where the kj are the exact eigenvalues or resonances of k, and rrij are their multiplicities. Eq. ( pi] ) gives the h expansion 

of (^) in terms of periodic orbit contributions. The first-order amplitudes A p °) of the expansion can be determined 
by adjusting the exact response function ( |l4| ) to the form of the semiclassical approximation 

g OBC (k)^Y, A Po eikSp ° ( 15 ) 

po 



by harmonic inversion [|25|. Note that for chaotic billiards the leading-order amplitudes A p ) as well as the higher- 
order amplitudes Apt) are independent of the wave number A;. In a direct harmonic analysis of the quantum spectrum 
only the zeroth-order term of the expansion ( pd| ) fulfills the ansatz for the harmonic inversion procedure. The higher- 
order terms act as a kind of weak noise which is separated from the zeroth-order "signal" by the harmonic inversion 
procedure. The direct harmonic analysis of the quantum signal will therefore yield exactly the lowest-order amplitudes 

The n th -order amplitudes ApV can be determined if the exact eigenvalues kj and their (n— l) st -order approximations 
kj tn -i are given. The (n — l) st -order approximation to the response function can be written in the form ([lj), with kj 
replaced with the approximation kj jH -i. One can then calculate the difference between the exact quantum mechanical 
and the (n — l) st -order response function and compare it with the expression resulting from the expansion (O), 



n— l oo oo 

5 qm (fc) - E 9m(k) = £ g m {k) = £ — £ A$ >e ia - fc . (16) 

m— m=n m—n po 

The leading-order terms in (|l^ ) are ~ k~ n , i.e., multiplication by k n yields 



9 qm 



= E4o )eiSpofe +°(T) • ( 17 ) 



po 



The right-hand side of (jT?]) now has assumed a form which is again suited to the harmonic inversion procedure. More 
precisely, the harmonic inversion of the weighted difference spectrum ( |l7[ ) will yield the periods s po and the n th -order 

amplitudes Apl) of the h expansion (HT 
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Fig. 1. Quantum mechanical density of states p(k) = ( — l/V) Im <?(fc) of the three-disk system with disk separation d — 6 (Ai 
subspace) as a function of real values of the wave number k. Only resonances of the four leading bands have been included. 
[Data courtesy of A. Wirzba.] 



4.2 Application to the three-disk scattering system 



We now apply the procedure to the three-disk system with disk separation d — 6. As a first step we perform a 
harmonic analysis of the exact quantum spectrum to obtain the leading-order (n = 0) periodic orbit contributions to 
the density of states. The exact quantum values for d = 6 were taken from Wirzba ^UgmBj. Figure |l| shows the 
quantum mechanical density of states p{k) = (— l/7r) Img(fc) for real values of the wave number k resulting from the 
four leading bands of the A\ subspace. Note that this set of resonances is of course not complete as the subleading 
bands with large negative imaginary part are not included. 

The spectrum in Figure [l] served as the signal for the harmonic inversion procedure. The results of the analysis 
turned out to be more accurate if the lowest part of the signal, determined by the "most quantum" resonances with 
very small real part, is cut off. We analyzed the spectrum in the range Rcfc £ [50,250] to obtain the periodic orbit 
contributions in two different intervals of the scaled action. The results are presented in Figure g. The solid lines 
give the sizes of the imaginary parts of the semiclassical amplitudes A^o calculated directly from Eq. ( [l2| ) using the 
classical periodic orbit data, as a function of the scaled action of the orbits. The crosses show the amplitudes resulting 
from the harmonic inversion of the quantum spectrum. Note the different scales of the two plots. The results of the 
harmonic inversion arc seen to be in excellent agreement with those from the classical calculations of the amplitudes 

^4po'' entering into Gutzwiller's trace formula, clearly confirming the validity of the latter. 

In a second step, we now determine the next-to-leading-order h corrections to the Gutzwiller trace formula for the 
three-disk system by the harmonic analysis of the difference spectrum between the exact quantum resonances and the 
semiclassical resonances of the A\ subspace. The semiclassical resonances for disk separation d — 6 had been calculated 
2(|^7j from a cycle expansion of the Gutzwiller- Voros zeta function. [Since the Gutzwiller- Voros zeta 
function is directly related to the Gutzwiller trace formula without further approximations, the semiclassical resonances 
resulting from both expressions will be the same.] The weighted difference spectrum is shown in Figure ||[ Note that 
due to the limited radius of convergence of the cycle expansion only resonances with Im k > —0.8 were available and 
could be included in the signal. 

In Figure I, the crosses designate the results for the first-order amplitudes A^o obtained from the harmonic inversion 
of the difference spectrum shown in Fig. 0, which was analyzed in the region Re k £ [100, 250]. For comparison, we also 

determined the first-order amplitudes A$ for each orbit following the method of Vattay and Rosenqvist described 
above (see Eq. ([l5|)). These results are represented in Figure || by solid lines. 



by Wirzba |21 
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Fig. 2. Imaginary parts of the amplitudes of the leading-order (n — 0) periodic orbit contributions to the density of states of 
the three-disk system with disk separation d = 6 as a function of the scaled actions of the symmetry reduced orbits. Solid lines: 



semiclassical amplitudes ^Ipo' versus scaled actions of the symmetry reduced orbits, calculated directly from classical mechanics. 
Crosses: results from the harmonic inversion of the exact quantum spectrum (Ai subspace). 



For almost all orbits, the harmonic inversion results for _4p are seen to be in perfect agreement with the amplitudes 
calculated by the method of Refs. [||, |[[l(]] . There is, however, one exception, namely the distinct discrepancy for the 
orbit with symbolic code '1' (scaled action s ss 4.267949). The deviation is systematic and appears in the same way 
if the parameters of the harmonic inversion procedure (such as signal length etc.) are varied. This point still needs 
further clarification. A possible explanation for the discrepancy may lie in the fact that the set of resonances from 
which the signal was constructed was not complete, since only resonances near the real axis could be included. However, 
this does not explain why only one orbit is strongly affected. On the other hand, the error might also be due to the 
theory of Refs. |p|j9[|ic|], or its application to the three-disk system. In fact, the V orbit is the orbit with the largest 
contributions from terms of higher order in I to the sum in (|13|) . The contributions from the different I terms and the 
converged sum over I of the five shortest orbits are given in Table [lj For comparison, the last column of Table [j] shows 
the corresponding values following from the amplitudes of the harmonic inversion results. The '1' orbit exhibits the 
largest deviation between the I = contribution and the converged sum over I, followed by the '0' orbit. For orbits 
with a symbol length of 2 or longer, the contributions of higher I terms to the amplitude A^o are already so small 
(due to the large absolute value of the stability eigenvalue A) that it is impossible to decide whether or not there 
is a discrepancy between these terms and the harmonic inversion results. (Note that this is also true for the period 
doubling of the orbits '0' and '1' in Table]!].) However, the harmonic inversion results for the '0' orbit, which also shows 
a relatively large contribution from the I = 1 term, are in agreement with the theory. Again, it cannot be explained 
why only the '1' orbit is affected (although in this case the reasons might lie in the special geometrical properties of 
the '0' orbit). 
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Fig. 3. Three-disk system with disk separation d = 6: Weighted difference spectrum k Ap(k) — k(p qm (k) — p sc (k)) between the 
quantum mechanical and the semiclassical density of states (Ai subspace) as a function of the wave number k. 

Table 1. Correction terms Cj (in units of the momentum hk) and their contributions /X 1 to the first-order h amplitude 
( p^[ ) for the five shortest periodic orbits of the three-disk system with d = 6. The values are compared with the results obtained 
by harmonic inversion (hi). The orbits are characterized by their symbolic code; their scaled action a and expanding stability 
eigenvalue A are also given. Note that the maximum correction to the I = contribution occurs for the orbit T', and is given 
by the I = 1 term. 





' Cl ~ 


oo „(1) 

\ ' °i 




"-3 ^(1) 


hi 


'0' 

s = 4.000000 
A = 9.898979 


0.625000 0.625000 

1 1.125000 0.113648 

2 -2.750000 -0.028064 

3 -14.750000 -0.015206 


0.690360 


0.6934 


T' 

s = 4.267949 
A = -11.77146 


1.124315 1.124315 

1 3.661620 -0.311059 

2 4.383308 0.031633 

3 1.162291 -0.000713 


0.843867 


1.055 


2x'0' 
s = 8.000000 
A = 97.98979 


1.250000 1.250000 

1 2.250000 0.022962 

2 -5.500000 -0.000573 

3 -29.500000 -0.000031 


1.272357 


1.259 


'01' 
s = 8.316529 
A = -124.0948 


2.039795 2.039795 

1 6.278740 -0.050596 

2 5.881196 0.000382 

3 -4.066328 0.000002 


1.989582 


2.019 


2xT' 
s = 8.535898 
A = 138.5672 


2.248630 2.248630 

1 7.323240 0.052850 

2 8.766615 0.000457 

3 2.324582 0.000001 


2.301937 


2.270 



8 



Kirsten Weibert et al.: Higher-order h corrections in the semiclassical quantization of chaotic billiards 



0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 



(a) 


X 


X 


x 

X 

X 


1 

X 



10 



12 



14 




Fig. 4. First-order h corrections to the trace formula of the three-disk system with disk separation d = 6 as a function of the 
scaled actions of the symmetry reduced orbits. Solid lines: first-order amplitudes following from a direct evaluation of Eq. (|l3|). 
Crosses: results from the harmonic inversion of the difference spectrum between exact quantum resonances and semiclassical 
cycle expansion values (Ai subspace). 



Concluding our discussion of Figure |j and Table we notice that the harmonic inversion results indeed confirm 
the validity of the I = approximation to the formula (||) for orbits with large stability eigenvalues. On the other 
hand, the results demonstrate that the theory of higher-order h corrections to the Gutzwiller formula still contains 
unanswered questions, and further investigations are necessary. 



5 Corrections to the semiclassical eigenvalues 



We now turn to the problem of obtaining corrections to the semiclassical eigenvalues of chaotic systems from the h 
expansion ( |ll|) of the periodic orbit sum. A general procedure for including higher-order h corrections in the harmonic 
inversion scheme was developed in Refs. ]lj,|l5],[ll|, where it was applied to the circle billiard as an example of an 
integrable system. In the following, we briefly recapitulate the main ideas of the procedure and then apply the technique 
to the open three-disk system. 



5.1 Theory 



For periodic orbit quantization, usually only the zeroth-order contributions ApJ to the expanded response function 
(|ll|), corresponding to the Gutzwiller formula (or, for integrable systems, the Berry- Tabor formula), are considered. In 
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the harmonic inversion scheme for semiclassical quantization |l^,[l3|,[L4|], semiclassical approximations to the eigenvalues 
or resonances are determined by adjusting the Fourier transform of the principal periodic orbit sum 

C (s)=Y,^S( S -s po ) (18) 

po 

to the functional form of the corresponding exact quantum expression (i.e., the Fourier transform of the exact response 
function (|l|)) 

C qm (s) = -iJ2 m J e ~ ifc,S , (19) 
j 

with kj the eigenvalues or resonances and rrij their multiplicities. 

Since for n > 1 the asymptotic expansion (|l"l| ) of the semiclassical response function suffers from the singularities 
at k = 0, higher-order h terms cannot directly be included in the harmonic inversion scheme. Instead, the correction 
terms to the semiclassical eigenvalues can be calculated separately, order by order. We assume that the (n — l) st -order 
ft, approximations fcj.„_i to the exact eigenvalues have already been obtained and the n th -order approximations kj iTl 
are to be calculated. In terms of these approximations to the eigenvalues, the difference between the two subsequent 
approximations to the quantum mechanical response function reads 



*.(*) = £ U 



kj,n + 10 k — fcj,n-l + iO 



with kj_ n = \{kj^ n + fcj, n _i) and Akj^ n = kj >n — fej-, n _i. Integration of (20) and multiplication by k n yields 

g n ( k ) = fc» j 9n (k)dk = y: T-f Ak +z ■ (21) 

3 3,n 



The periodic orbit approximation to ( |2l| ) is obtained from the term g n (k) in the periodic orbit sum ( |ll[ ) by integration 
and multiplication by k n , yielding 

g n (k) = -iJ2— Ai n M k »>° + o(l) . (22) 

po S P° 

One can now Fourier transform both ( pl| ) and (p^), and obtains (n > 1) 

Cn{s) = — / g n (k)e- [ks dk 



2tt _ 

= i£m i (fc i )"Z\fc J ,„e- i ^ s (23) 

i 

= -iE— 4o ) *(*-*po). (24) 

po S P° 

Equations (p3| ) and (^i|) imply that the H expansion of the semiclassical eigenvalues can be obtained, order by order, 
by adjusting the periodic orbit signal (53) to the functional form of ( p3| ) by harmonic inversion (h.i.). The frequencies 
kj of the periodic orbit signal (|2J) are the semiclassical eigenvalues or resonances, averaged over different orders of h. 
Note that the accuracy of these values does not necessarily increase with increasing order n. The corrections Akj^ n to 
the eigenvalues are not obtained from the frequencies, but from the amplitudes, rrij(kj) n Akj_ n , of the periodic orbit 
signal. 



5.2 Application to the three-disk scattering system 

We have applied the above procedure to the open three-disk scatterer with disk separations d — 6 and d = 2.5. For 
both cases, we first calculated the zeroth-order H approximations to the resonances and then determined the first-order 
h corrections to the semiclassical results following the scheme outlined above. 
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For disk separation d = 6, we used the periodic orbits up to length s max = 56 to calculate the resonances in the 
region < Refc < 250. In the first-order amplitudes (|l3|), only the leading-order term Z = of the sum was included. 
(The terms with I > 1 contribute significantly only to the '1' orbit and thus basically do not effect the semiclassical 
resonances.) The results for the first-order corrections Aki were added to the semiclassical results to obtain the 
first-order approximations k\ to the resonances. Table ^ shows part of the results in the regions Re k 6 [0, 12] and 
Refc G [150, 155] with Imfc > —0.5. For comparison, the exact quantum resonances k ex from Refs. (2l],^6|,^7j are also 
given. We note that the semiclassical resonances kg obtained by harmonic inversion agree with the semiclassical cycle 
expansion values from Refs. [|l],|2(||27j to all digits given. Figure [| compares the semiclassical errors of the zeroth-order 
(crosses) and first-order (squares) approximations as a function of the real part of the resonances. The deviations of 
the real and imaginary parts from the exact quantum values are shown separately. Only resonances with an imaginary 
part Irak > — 0.5 were included in the plot. 

The results presented in Figure [s] show that by including the first-order corrections a significant improvement in 
the accuracy of the real parts of the resonances is achieved. This is evident from Figure g in spite of the fact that 
a one-to-one correspondence between the zeroth and first-order values plotted is difficult to establish with the naked 
eye. For most resonances, the real part of the first-order approximation lies between two and five orders of magnitude 
closer to the exact quantum values than the zeroth-order approximation. Only for the "most quantum" resonances, 
with very low real parts, the improvement is rather small. The reason for this lies in the nature of the semiclassical 
approximation as an approximation itself: in order to improve these values, second or higher-order terms of the H 
expansion must be considered. 

The accuracy of the imaginary parts of the semiclassical resonances is less significantly increased by the first-order 
corrections than that of the real parts. For some resonances, the zeroth-order approximation lies even closer to the 
exact quantum values than the first-order approximation. This was also observed in Refs. |6|[l(|, where the first-order H 
corrections to the resonances was calculated using the cycle expansion technique. As discussed in |||l0|, the first-order 
corrections to the periodic orbit sum mainly improve the real part of the resonances, while the imaginary part can be 
expected to be improved by second-order H corrections. 

A similar behaviour can be found for disk separation d = 2.5. Here, we calculated the semiclassical resonances 
and their first-order h corrections in the range < Re/c < 90 and —0.82 < Imk < from the periodic orbits up to 
length s max = 12. In the first-order amplitudes (|l3|), again only the I — term was included. Table |^ compares part 
of the results for the first-order approximations to the resonances with the zeroth-order approximations and the exact 
quantum values. The zeroth order resonances ko obtained by harmonic inversion agree with the cycle expansion values 
calculated by Wirzba [^l,^6|,[27j] (not shown) to at least four significant digits. 

Again, we determined the semiclassical error of the first-order approximations to the resonances in comparison to 
that of the zeroth-order approximation. The results are presented in Figure ||. The general behaviour of the values is 
similar to that in the case d = 6 discussed above, although the improvement of the accuracy achieved by the first-order 
corrections is not as spectacular as for d — 6. The reason for this may partly lie in the error induced by the harmonic 
inversion method, which for d — 2.5 is larger already in the zeroth-order approximation than for d = 6. The results 
could in principle be improved by extending the signal to longer orbits. On the other hand, in the part of the spectrum 
considered, second and higher-order h corrections may be more important than in the case d = 6. However, it is evident 
from Figure ^| that, apart from the resonances with very small real parts, the semiclassical error of the real parts of 
the resonances could still be reduced, even for the small disk separation of d = 2.5, by the first-order h corrections by 
one or two orders of magnitude. 



6 Conclusions 

In this paper we have demonstrated the power of the harmonic inversion technique in explicitly determining higher- 
order h corrections to the Gutzwiller trace formula and to the semiclassical eigenvalues of a completely chaotic system, 
namely the three-disk scattering system. The method has been used in two directions: (1) for the harmonic analysis of 
the exact quantum spectrum, (2) for the direct calculation of higher-order corrections to the semiclassical eigenvalues 
from classical periodic orbit data. 

The harmonic analysis of the exact quantum spectrum of the three-disk system with the "standard" literature 
disk separation of d = 6 first yielded the zeroth-order semiclassical amplitudes of the periodic orbit sum (i.e., the 
amplitudes entering the Gutzwiller formula) , which were found to be in perfect agreement with the Gutzwiller ampli- 
tudes calculated directly from classical periodic orbit data. Next, from the exact quantum resonances and their zeroth 
order approximations, we were able to compute the first-order amplitudes applying harmonic inversion to Eq. (|l7|). 
We could verify the correctness of the values obtained in this way by comparing with the results of an alternative 

] for calculating first-order corrections to the Gutzwiller formula in chaotic billiards, which 
we implemented for the three-disk system. The results turned out to be in very good agreement (on the order of 1.5 
per cent, or better), with one notable exception, namely the '1' orbit, for which a distinct discrepancy (on the order 
of 20 per cent) persisted. We have discussed possible origins of the discrepancy although an ultimate reason could not 
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Kirsten Weibert et al.: Higher-order h corrections in the semiclassical quantization of chaotic billiards 



11 



Table 2. Zeroth (fco) and first (fci) -order approximations to the complex eigenvalues of the resonances of the three-disk system 
with disk separation d — 6 (Ai subspace), obtained by harmonic inversion of a signal of length s max = 56. For comparison, the 
exact quantum values fc cx are given (taken from Refs. [ pH|2^ ,p7|). Only resonances of the four leading bands with Imfc > —0.5 
are included. 
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Table 3. Zeroth and first-order approximations to the complex eigenvalues of the resonances of the three-disk system with disk 
separation d = 2.5 (Ai subspace), obtained from a signal of length s max = 12. The notations are the same as in Table ^. The 
table contains the resonances in the region 1 < Refc < 90 and —0.82 < Imfc < 0. 
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Fig. 5. The semiclassical errors of the zeroth (+) and first (□) -order approximations to the complex eigenvalues of the 
resonances of the three-disk system with d = 6, plotted as a function of the real part of the resonances. Only resonances with 
imaginary parts Imfc > —0.5 arc included. 
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Fig. 6. As Fig. ^, but for disk separation d — 2.5. Only resonances with imaginary parts Imk > —0.82 are included. 



be identified. Therefore, in spite of the very good agreement of the results in all other cases, we have to conclude that 
the theory of h corrections to the Gutzwiller formula still contains unanswered questions. [We note that in fact for 
integrable systems there does not yet exist a general theory for higher-order h corrections to the Berry- Tabor formula 
at all.] 

In the direct calculation of higher-order corrections to semiclassical eigenvalues from classical periodic orbit data, 
we first evaluated the first-order correction amplitudes to the Gutzwiller formula as given by the theory of Vattay and 
Rosenqvist, and then, by harmonic inversion of Eq. J24|), determined the first-order H corrections to the semiclassical 
(complex) eigenvalues of resonances of the three-disk scattering system with disk separations d — 6 and d = 2.5. For 
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both distances, the semiclassical error, as compared to the exact quantum values, of the zeroth-order results (obtained 
from the Gutzwiller formula by harmonic inversion) for the real parts of the resonances could be significantly reduced 
by including the first-order h corrections: the accuracy was increased by two to five orders of magnitude for d = 6, 
and still by one to two orders of magnitude for d — 2.5. Only for the "most quantum" resonances, with very small 
real parts, the increase in accuracy was found to be rather modest. It turned out that the accuracy of the imaginary 
parts of the semiclassical eigenvalues of resonances was less significantly increased by the first-order corrections; here, 
second-order corrections would have to be considered. 

Although in our calculations we have used literature values for the semiclassical resonances and for the zeroth and 
first-order amplitudes, we could have performed, in principle, all calculations knowing the exact quantum resonances 
only. The analysis of the exact quantum spectrum yields the semiclassical amplitudes, which in turn can be used to 
calculate the semiclassical resonances. Then, by an analysis of the difference spectrum between semiclassical and exact 
resonances, the first order amplitudes can be determined, which again can be used to obtain the first-order corrections 
to the resonances. Although we have concentrated in our examples on obtaining first-order h corrections, it is evident 
from our discussion that the next-order corrections could be obtained iteratively in an analogous manner by repeated 
application of Eqs. (|l7|), and (|23|), (p4|), respectively. 

In summary, we have demonstrated that harmonic inversion - as a means for circumventing the convergence 
problems of semiclassical trace formulae - is indeed a very efficient and universal tool, not only for semiclassical 
quantization, but also for the explicit calculation of higher-order H corrections to the semiclassical eigenvalues or 
resonances even in chaotic systems. Moreover, harmonic inversion does not rely on specific assumptions for the systems 
under consideration, and therefore an application of the methods presented in this paper to other chaotic systems will 
be promising and worthwhile. 

We are grateful to G. Vattay for sending us the code for the calculation of the first-order corrections to the Gutzwiller formula 
for chaotic billiards. We thank A. Wirzba for communicating to us quantum mechanical and semiclassical resonances of the 
three-disk system. This work was supported by Deutsche Forschungsgemeinschaft and Deutscher Akademischer Austauschdienst. 
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